Nonequilibrium phase transitions and finite size scaling 
in weighted scale-free networks 
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We consider nonequilibrium phase transitions in weighted scale-free networks, in which highly 
connected nodes, which are created earlier in time are partially immunized. For epidemic spreading 
we solve the dynamical mean-field equations and discuss finite-size scaling theory. The theoretical 
predictions are confronted with the results of large scale Monte Carlo simulations on the weighted 
Barabasi- Albert network. Local scaling exponents are found different at a typical site and at a node 
with very large connectivity. 
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I. INTRODUCTION 

Complex networks, which have a more complicated 
topology than periodic lattices have been observed in a 
large class of systems in different fields of science, tech- 
nics, transport, social and political life, etc, see Refs i 1 ' 2 ' 3 
for recent reviews. The structure of complex networks is 
described by graphs^, in which the nodes represent the 
agents and the edges the possible interactions. Generally 
there are connections between remote sites, too, which 
is known as the small world effect. Another feature of 
many real networks is the non-democratic way of the dis- 
tribution of the links: there are sites which are much 
more connected than the average and the distribution of 
the number of edges, P(k), has a power-law tail: 



P D {k) ~ Ak- 



k > 1 , 



(1) 



thus the edge distribution is scale free. In real networks 
the degree exponent is generally: 2 < 7 < 3 and as shown 
by Barabasi and Albert scale-free networks are usually 
the results of growth and preferential attachment. 

Since the agents of a network interact in one or an- 
other way it is natural to ask about the cooperative be- 
havior of the system. In particular if there exist same 
kind of (thermodynamical) phases and if there are sin- 
gular points as the strength of the interaction or other 
suitable parameter (such as the strength of disordering 
field, temperature, etc.) are varied. In this respect static 
models&SiiSiii (Ising, Potts models, etc.), as well as non- 
equilibrium processe a 12 i 13 i 14 (percolation, spread of epi- 
demics, etc.) are investigated. Generally non-weighted 
networks are considered, in which the strength of inter- 
action at each bond is constant. Due to long-range in- 
teractions conventional mean-field behavior is expected 
to hold, at least if the network is sufficiently weakly con- 
nected, i.e. 7 > 7„. For the Ising model 7 U = 5, whereas 
for the percolation and epidemic spreading j u — 4. At 
7 = 7m there are logarithmic corrections to the mean- 
field singularities, whereas for 7^ > 7 > "f c we arrive to 
the unconventional mean-field regime in which the criti- 
cal exponents are 7 dependent. The effect of topology of 
scale- free networks becomes dramatic for 7 < 7 C , when 



the average of k 2 , (k 2 ), as well as the strength of the aver- 
age interaction becomes divergent. Consequently for any 
finite value of the interaction scale-free networks are in 
the ordered state, c.f. there is no threshold value of epi- 
demic spreading. Since 7 C = 3, in realistic networks with 
homogeneous interactions always this type of phenomena 
should occur. 

Recently, much attention has been paid on weighted 
networks, in which the interactions are not homogeneous. 
Generally the strength of interactions of highly connected 
sites are comparatively weaker than the average, which 
can be explained by technical or geographical limitations. 
To model epidemic spreading one should keep in mind 
that sites with a large coordination number are generally 
earlier connected to the network and in the long period 
of existence they have larger chance to be immunized. 

An interesting class of degraded networks has been 
introduced recently by Giuraniuc et aU^ in which the 
strength of interaction in a link between sites i and j is 
rescaled as: 



(kikj 



(2) 



where fcj and kj are the connectivities in the given sites. 
The properties of this type of network in equilibrium crit- 
ical phenomena, in particular for the Ising model have 
been studied in detail in Re£i£. Most interestingly the 
equilibrium critical behavior is found to depend on the 
effective degree exponent, 



7' = (7 - /i)/(l - M) 



(3) 



thus topology and interaction seem to be converted. One 
important aspect of the degraded network in Eq. J2J that 
phase transition in realistic networks with 7 < 3 is also 
possible, if the degradation exponent, /1, is sufficiently 
large. Therefore theoretical predictions about critical 
singularities can be confronted with the results of nu- 
merical calculations. 

In this paper we study nonequilibrium phase transi- 
tions in weighted networks. Our aim with these in- 
vestigations is twofold. First, we want to check if the 
simple reparametrization rule in Eq.(|3J) stays valid for 
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nonequilibrium phase transitions, too. For this purpose 
we make dynamical mean-field calculations and perform 
large scale Monte Carlo simulations. Our second aim is 
to study the form of finite-size scaling in nonequilibrium 
phase transitions in weighted scale-free networks. To an- 
alyze our numerical results we use recent field-theoretical 
calculations 7 in which finite-size scaling in Euclidean lat- 
tices above the upper critical dimension has been stud- 
ied. In the conventional mean-field regime of scale-free 
networks analogous scaling relations are expected to ap- 
ply- 

The structure of the paper is the following. The dy- 
namical mean-field solution of the problem is presented 
in Secl^ whereas finite-size scaling theory is shown in 
Sec lIIII Results of Monte Carlo simulations of the con- 
tact process on weighted Barabasi- Albert networks are 
presented in Sec II VI and discussed in SeclVl 



the dynamical mean- field equations in Eq.Q) are given 
by: 



dt 



= ~ Xk l M (! - Pi)P - Pi 



(7) 



with A = A < k 1 -* 1 > I < k >< k^ 11 > 2 . In the sta- 
tionary state, dpi/dt — 0, the local densities are given 
by: 



Afc. 1 " 



Pi 



1 + Xk}-* 1 



P 



(8) 



i.e. they are proportional to k\ M . Putting pi from Eq. JHJl 
into Eq.© we obtain an equation for p: 



< k 1 -" >= A 



P D {k) 



fc2(l-M) 

1 + xk}-^, 



-dk 



(9) 



II. DYNAMICAL MEAN-FIELD SOLUTION 

In the calculation we consider the contact process^, 
which is the prototype of a non-equilibrium phase tran- 
sition in the directed percolation universality classi 7 -. In 
this process site i of the network is either vacant (0) or 
occupied by at most one particle (A). The dynamics of 
the model is given by a continuous time Markov process 
and is therefore defined in terms of transition rates. The 
reactions in the system are of two types: a) branching 
in which particles are created at empty sites (provided 
one of its neighbors, j, is occupied) occurs with rate Xij, 
b) death of particles with a rate, p. This latter rate we 
set [i = l. In this section we solve the problem in the 
mean- field approximation, which is expected to be exact, 
due to long-range interactions in the system. 

We start with the set of equations for the time deriva- 
tive of the mean- value of the density, pi, at site, i = 
1,2 N: 



dp 



=^2 X i,j( 1 - Pi) Pi 



Pi 



(4) 



and correlations in the densities at different sites are 
omitted. In the next step in the spirit of the mean-field 
approach we replace the interactions: 



_ kjkj (kjkj) *- 



(5) 



i.e. there is an interaction between each site the (mean) 
value of which is proportional to the probability of the 
existence of that bond. Now in terms of an average den- 
sity: 



Pi 



51 j k 3 



(G) 



where summation over i is replaced by an integration 
over the degree distribution, Pr>(k), and in the ther- 
modynamic limit the upper limit of the integration is 
kmax — * oo. The solution of EqJHJ in the vicinity of the 
transition point, p< 1, depends on the large- A; limit of 
the degree distribution in Eq. and is given in terms of 
the integration variable, k' = fc 1_Al , as: 

< k' >= XA ( jfe'-T' — K — dfc' = Q{p,i), p<l 



1 + Xk'p 



(10) 

where 7' is defined in Eq.©. Note, that the functional 
form of the equation in (| 1 0|> is identical to that for stan- 
dard scale-free networks, just with an effective degree 
exponent, 7'. Thus the solution in Ref^ can be applied 
and in this way we have obtained an extension of the 
results in Ref^ for nonequilibrium phase transitions. 

To analyze the solution of Eq. (|1U|I we apply the method 
in Ref. — , which is somewhat different from the original 
method in Refiii 

• 7' > 4 

For small p, Q(p, 7') in Eq. (fTU|l can be expanded 
in a Taylor series at least up to a term with ~ 
p . Consequently there is a finite transition point, 
A c =< k' > j A < k' 2 >, and the density in the 
vicinity of the transition point behaves as: p(X) ~ 
(A c — A). This is the conventional mean-field regime. 
At the borderline case, 7' = 4, there are logarithmic 
corrections to the mean- field singularities. 

• 3 < 7' < 4 

For small p only the linear term in the Taylor ex- 
pansion of Q{p, 7') exists. The p-dependence of the 
next term, ai (p), is singular and given by: 



02 (p) = -X 2 pA 



k'-i 



Xk'p 



-dk' 



(11) 
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The p-dependence can be estimated by noting that 
for a small, but finite p there is a cut-off value, 
k' ~ 1/p, so that 

a 2 (p) ~ -\ 2 pA J k'-^'k^dk' ~ p 7 '" 3 . (12) 

Consequently the density at the transition point 
behaves anomalously, 

p~(\ c -\f, /3 =l/( 7 ' -3). (13) 

This is the unconventional mean- field region. 

• i < 3 

In this case Q(p, 7) is divergent for small p. Its 
behavior can be estimated as in Ea. 1)12(1 leading to 
Q{p>l') ~ A 7 ~ 2 p 7 ~ 3 . Consequently the system 
for any non-zero value of A is in the active phase. 
As A goes to zero the density vanishes as: 

p ~ A (7'-2)/(3- 7 ') _ (w) 

Here at the border, 7' = 3, the system is still in the 
active phase, but the density is related to a small 
A as: |ln(pA)| ~ 1/A. 

Before we confront these analytical predictions with 
the results of numerical simulations we discuss the form 
of finite-size scaling in scale-free networks. 

III. FINITE-SIZE SCALING 

In a numerical calculation, such as in Monte Carlo 
(MC) simulations, one generally considers systems of fi- 
nite extent and the properties of the critical singulari- 
ties are often deduced via finite-size scaling. It is known 
in the phenomenological theory of equilibrium critical 
phenomena that due to the finite size of the system, 
L, critical singularities are rounded and their position 
is shifted^. As it is elaborated for Euclidean lattices 
finite-size scaling theory has different forms below and 
above the upper critical dimension, d c . For d < d c in 
the scaling regime the singularities are expected to de- 
pend on the ratio, L/£, where £ is the spatial correla- 
tion length in the infinite systemic. On the other hand 
for d > d Cl when mean-field theory provides exact val- 
ues of the critical exponents, finite-size scaling theory 
involves dangerous irrelevant scaling variables^, which 
results in the breakdown of hyperscaling relations. For 
equilibrium critical phenomena predictions of finite-size 
scaling theory^! above d c are checked numerically, but 
the agreement is still not satisfactory^. 

For non-equilibrium critical phenomena finite-size scal- 
ing above d c has been studied only very recently^ and 
here we recapitulate the main findings of the analysis. 
For directed percolation, which represents a broad class 
of universality^, dangerous irrelevant scaling variables 



are identified in the fixed point. As a consequence scal- 
ing of the order-parameter is anomalous: 

p = L-^»* piSL 1 ^' ,hL A ^*) , (15) 

Here, S, is the reduced control parameter, with the nota- 
tions of Sec. [H](5 = (A — A c )/A c and h is the strength of an 
ordering field. The critical exponents, f3 = 1 and A = 2, 
are the same as in conventional mean-field theory. The 
finite-size scaling exponent is given by, v* — 2/d, and 
thus depends on the spatial dimension, d. Note, that be- 
low d c = 4 it is the correlation length exponent, v, which 
enters into the scaling expression in Ea. ll5fl . but above 
d c , due to dangerous irrelevant scaling variables it should 
be replaced by v* . At the critical point, 8 — 0, the scal- 
ing function, p(0, x), has been analytically calculated and 
checked by numerical calculations. 

In the following we translate the previous results for 
complex network, in which finite-size scaling is naturally 
related to the volume of the network, which is given by 
the number of sites, N. In the conventional mean-field 
regime with the correspondence, N <-> L d , we arrive from 
Eci. (|15|) to the finite-size scaling prediction: 

Pt yP = N-^ 2 p typ (SN^ 2 ,hN A / 2 ) , (16) 

which is expected to hold for a typical site, i.e. with 
a coordination number, k ~ (fc). On the other hand 
for the maximally connected site with k max ~ iV 1 ^ 7-1 ) 
according to Eq. JSJ the finite-size scaling form is modified 
by: 

p max = N-V^-M^p max {SN x l 2 ,hN A / 2 ) . (17) 

Since in the derivation of the relation in Eq. I|15|l the ac- 
tual value of (3 has not been used, we conjecture that the 
results in Eas. (|16|> and (|17(l remain valid in the uncon- 
ventional mean-field region, too. 

IV. MONTE CARLO SIMULATION 

In the actual calculation we considered the contact pro- 
cess on the Barabasi- Albert scale- free network 6 , which 
has a degree exponent, 7 = 3, and used a degradation ex- 
ponent, p — 1/2. Consequently from Eq.© the effective 
degree exponent is 7' = 5, thus conventional mean-field 
behavior is expected to hold. (We note that the same 
system is used to study the equilibrium phase transition 
of the Ising model in Refii^.) Networks of sites up to 
N = 1024 are generated by starting with mo = 1 node 
and having an average degree: (k) = 2. Results are aver- 
aged over typically 10000 independent realizations of the 
networks. 

In the calculation we started with a single particle at 
site i, (which was either a typical site or the maximally 
connected site) and let the process evolve until a station- 
ary state is reached in which averages become time inde- 
pendent. In particular we monitored the average value of 
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FIG. 1: Variation of the order parameter at a typical site as 
a function of the creation rate, A. Inset: Enlargement in the 
critical region. 




FIG. 2: As in Fig^fc* the maximally connected site. 



the occupation number, pi, as introduced in mean- field 
theory in Eq.(@J, and the fraction of occupied sites, m,-, 
(order parameter) as a function of A, whereas p was set 
to be unity. In the stationary state and in the vicinity 
of the transition point pi and m s are expected to be pro- 
portional with each other and characterize the order in 
the systemSi. 

The A dependence of the order parameter is shown in 
Figffl f° r a typical site and in Fig [3 for the maximally 
connected site. Evidently there is a phase transition in 
the system in the thermodynamic limit, which is rounded 
by finite size effects as shown in the insets of Figs, ^and 

121 

To locate the phase transition point we form the ra- 
tios: r(N) = m(N)/m(N/2), for different finite sizes. As 
shown in Fig[3|r(iV) tends to zero in the inactive phase, 
A < A c , and tends to a value of one in the active phase, 
A > A c . The curves for different N cross each other and 
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FIG. 3: The ratio r(N) = m(N) /m(N/2) a) for a typical site, 
b) for the maximally connected site. Note that the location 
of the crossing points, which defines A c is the same in the two 
cases, whereas the value at crossing, which is related to the 
finite-size exponent, x, through r(N,X c ) = 2~ x is different. 



the crossing point can be used to identify A c through ex- 
trapolation. Furthermore the value of the ratio at the 
critical point is given by: r(N, A c ) = 2~ x , where x is the 
finite-size scaling exponent, as given in scaling theory in 
Eqs.OU and |T7|). 

The transition point is found to be the same within 
the error of the calculation both for a typical site and 
for the maximally connected site and given by: A c = 
2.30(1). The finite-size scaling exponent, however, calcu- 
lated from r(N, A c ) is different in the two cases. For a 
typical site we estimate: Xt yp = 0.54(5), which should be 
compared with the field-theoretical prediction in Ea. (|16fl . 
which is x t y P = (3/2 — 1/2. In the maximally con- 
nected site the finite-size scaling exponent is measured as 
Xmax — 0.27(3), which again agrees well with the field- 
theoretical prediction in Ea. (|T7|) : x max — /3/2 — (1 — 
M )/(7-l) = l/4. 

Next, we consider correlations in the vicinity of the 
transition point and calculate the relation between the 
correlated volume, V, and the distance from the critical 
point, <5, which is expected to be in a power-law form, V ~ 
(5|~ w . According to field-theoretical results in Eas. i|lt)|) 
and i|17[l this exponent is uo — 2, both at a typical site 
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FIG. 4: Scaling collapse of the order-parameter near the 
transition point using the functional forms in Eg. 1161 for a 
typical site (a) and in Eg. 11711 for the maximally connected 
site (b). The finite-size scaling exponents, x typ and i miB are 
fixed by the previous analysis and the correlation exponent, 
lo, is used to have an optimal collapse of the data, see text. 



and at the maximally connected site. Now in the limiting 
case, V ~ N, the scaled order-parameter, rh = mN x , is 
expected to depend on the scaling combination, N X ^S, 
which is demonstrated in Fig. 0] both in the typical site 
(a) and in the maximally connected site (b). In both 
cases x and A c are fixed by the previous analysis and lo 
is obtained from the optimal scaling collapse, as Lo typ = 
2.05(10) and uj ma x = 2.00(5). Thus, once more we have 
a good agreement with the field-theoretical results. 

Finally, we turn to analyze dynamical scaling in the 
system. At the critical point, A = A c we have mea- 
sured the number of active sites, N a , as a function of 
time, t, which is shown in Fig|5]m a log-log plot, when 
the starting point is a typical site or the maximally con- 
nected site. As seen in FigO this relation is asymptot- 
ically given by N a ~ i a , where the effective value of 
the critical exponent, a, has a strong size dependence, 
in particular by starting with a typical site. By extrap- 
olation we obtained at yp = 0.98(5), which is compati- 
ble with the mean-field and finite-size scaling prediction, 
o-typ — 1, whereas starting from the maximally connected 
site we extrapolated a max = 0.57(2). Now, keeping in 
mind that the correlated volume can be expressed as, 

we obtain the relation 
t ~ VS with ( = (1 — x)/a. For a typical site our 
numerical result, Ctyp = 0.47(5) is compatible with the 
theoretical prediction of Ctyp = 0.5, whereas for the max- 
imally connected site we obtained Cmax = 1.28(5). From 
these results, using N ~ <5~ w , we obtain for the scal- 
ing behavior of the relaxation time, r, in the vicinity of 
the transition point, r ~ §- v } with v 1 - = (oj, so that 
v typ — 0.96(5) and v^ ax = 2.5(1). Note that once again 
at the typical site we are in complete agreement with the 



FIG. 5: Time dependence of the number of active sites start- 
ing from a single seed (a) a typical site, b) the maximally 
connected site) at the critical point in a log-log scale, Due to 
finite-size effects the curves are saturated for large time. The 
slope of the straight part of the curves defines the effective 
exponent, a e ff, which has strong finite-size dependence. We 
draw the staight lines for the largest size having a) a e ff = .8 
and b) a e ff = -5. 



mean- field and finite-size scaling result, vh p = 1. 



V. DISCUSSION 

We considered non-equilibrium phase trasnsitions in 
weighted scale-free networks, in which the creation rate 
of particles at given sites is rescaled with a power of the 
connectivity number. In this way non-equilibrium phase 
transitions arc realized even in realistic networks having 
a degree exponent, 7 < 3. Mean-field theory, which is 
generally belived to be exact in these lattices, is solved 
and the previously known three regimes of criticality 
(conventional and unconventional mean-field behavior, as 
well as only active phase) are identified. The theoreti- 
cal predictions in the conventional mean-field regime are 
confronted with the results of Monte-Carlo simulations 
of the contact process on the weighted Barabasi-Albert 
network. To analyze the simulation results we have 
applied and generalized recent field-theoretical results- 
about finite-size scaling of non-equilibrium phase tran- 
sitions above the upper critical dimension, i.e. in the 
mean-field regime. For a network the natural variable is 
the volume (mass) of the system which enters in a sim- 
ple way into the scaling combinations. We have obtained 
overall agreement with this finite-size scaling theory in 
which the critical exponents are simple rational numbers. 
We have also numerically demonstrated that at sites with 
very large connectivity there are new local scaling expo- 
nents, which differ from the values measured at a typical 
site. 
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